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Abstract: Multiple partonic interactions are widely used to simulate the hadronic 
final state in high energy hadronic collisions, and successfully describe many features 
of the data. It is important to make maximum use of the available physical con- 
straints on such models, particularly given the large extrapolation from current high 
energy data to LHC energies. In eikonal models, the rate of multiparton interactions 
is coupled to the energy dependence of the total cross section. Using a Monte Carlo 
implementation of such a model, we study the connection between the total cross 



section, the jet cross section, and the underlying event. By imposing internal consis- 
tency on the model and comparing to current data we constrain the allowed range of 
its parameters. We show that measurements of the total proton-proton cross-section 
at the LHC are likely to break this internal consistency, and thus to require an ex- 
tension of the model. Likely such extensions are that hard scatters probe a denser 
matter distribution inside the proton in impact parameter space than soft scatters, 
a conclusion also supported by Tevatron data on double-parton scattering, and/or 
that the basic parameters of the model are energy dependent. 
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1. Introduction 

Hadron-hadron collision events at high energies often contain high transverse energy 
jets, which in QCD arise from gluon or quark (generically, parton) scattering followed 
by QCD radiation and hadronization. This model is generally taken to be realistic 
above some minimum transverse momentum scale, p™ m . The contribution of these 
events to the total cross section rises with hadron-hadron centre-of-mass energy, s, 
since the minimum value of the x probed is given by 4(p™ in ) 2 /s, and the parton 
densities rise strongly for x < 10~ 2 or so [1,2]. 

One reason that this rising contribution to the cross section is of interest is that 
while perturbative QCD cannot predict total hadronic cross sections (since in many 
events no hard perturbative scale is present), rising hadronic cross sections are a 
feature also seen in successful non-perturbative approaches [3,4], the behaviour of 
which must presumably emerge from the QCD Lagrangian in some manner. Thus 
by comparing the behaviour of the hard contribution to the cross section to the 
behaviour expected from fits to the total cross section, it may be possible to learn 
something about the connection between these approaches and about hadronic cross 
sections in general. 
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The connection between the hard partonic cross section and the total cross sec- 
tion is not one-to-one, however. There are certainly hadronic scatters in which no 
hard jets are produced, and some non-perturbative scattering process must be added 
to the perturbative jet contribution to model the total cross section. In addition, 
at the high parton densities probed at recent, current and future colliders, simple 
assumptions lead to the conclusion that the probability of multiple partonic scatters 
in a single hadron-hadron collision is significant. In fact, Fig. [TJ shows that for pf 1111 
values below about 5 GeV, the total "hard" cross section calculated assuming one 
parton-parton scatter per proton-proton collision exceeds the total cross section as 
extrapolated using the non-perturbative fits, at LHC energies. This strongly implies 
that the average number of partonic scatters in an inelastic collision must be greater 
than one. 

Introducing the possibility of 
such multiparton interactions also 
seems to be required in order to 
describe the hadronic final state [5- 
7]. In general, softer additional 
scatters occurring in a high-p t 
event manifest themselves as ad- 
ditional particles and energy-flow, 
the so-called "underlying event" . 

In this paper we examine the 
predictions of the model that was 
discussed, for example in [8], and 
implemented in [9-12] including 
the possibility of soft scatters. We 
explore the consistency constraints 
that would be imposed by compar- 
ing a given value of the total cross 
section to the predicted jet cross 
section, and attempt to identify al- 
lowed regions of parameter space 
within which the model must lie if it is to be consistent with the measured cross 
section at the LHC. We also discuss ways in which energy dependencies in the pa- 
rameters could arise, and their impact upon these constraints. The studies are all 
carried out using the new implementation in Herwig++ [11,12]; however, they are 
also relevant to the fortran implementation Jimmy [9], if the same hard cross section 
is used. 
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Figure 1: The inclusive hard cross section for three 
different proton PDFs, compared to various extrap- 
olations of the non-perturbative fits to the total pp 
cross section at 14 TeV centre-of-mass energy. 
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2. Total and elastic cross section parameterizations 



Throughout this paper we will exploit the connection that can be established between 
the eikonal model of Refs. [9-11] and the total cross section. To give a reasonable 
range of expectations for the latter, we use the successful parameterization of Don- 
nachie and Landshoff [3,4]. We will use three different variations; 



1. The standard parameterization from [3] with the following behaviour at high 
energies: 

i 0.0808 



a tot ~ 21.7 mb • (-7^772) -> ^(U TeV) = 101.5 mb . (2.1) 



2. Using the same energy dependence but normalizing it to the measurement [13] 
by CDF: 

/ o \ 0.0808 

a tot ~ 24.36 mb • (-^) -> a tot (14 TeV) = 114.0 mb . (2.2) 

3. Using the most recent fit [4], which takes the contributions from both hard and 
soft Pomerons into account: 

/ o \ 0.0667 / s \ 0.452 

a tot ~ 24.22 mb • ( 5 ) + 0.0139 mb • I 5 ) 

V GeV 2 / V GeV 2 / (2.3) 

-> cr tot (14 TeV) = 164.4 mb . 

Other parameterizations and models for the total cross section exist [14,15], but 
their predictions for the total cross section at 14 TeV generally lie within the range 
covered by these three 1 . As will be seen, the range is wide, and early measurements 
of the total cross section at the LHC can be expected to have a big impact [18]. 

We will also find it useful to compare our model with the elastic slope parameter, 
B, defined in terms of the differential elastic scattering cross section, dcr/dt, as 



B = B(s,t = 0) 



(2.4) 



J t=o 



In the Donnachie-Landshoff parameterization, this is given by: 



B = 2a In — + B (2.5) 
so 



with a' = 0.25 GeV' 2 . Together with the CDF data [19], this implies 

B=(ln ^ — + (17 ±0.25)) GeV' 2 = (In ^ + (19 ± 0.25) ) GeV" 

V 1800 GeV V 7 V 14 TeV / 



■2 



(2.6) 



1 Thc most recent models of [16,17] predict (T to t(14 TeV) ~ 90 mb, which is 10 % below the 
smallest expectation we use. Since the difference this introduces is similar to the one between our 
first and second parameterization it can easily be estimated by the reader. 
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The most recent fit [4] has the same value for a' and hence B, while those of [14, 15] 
are a little higher: 20-22 GeV~ 2 . We therefore use the CDF value for the Tevatron 
energy and the range 19-22 GeV -2 to represent the range of possible measurements 
from the LHC. 



3. Eikonal model 

The scattering amplitude A(s, t) can be expressed as the Fourier transform of the 
elastic scattering amplitude a(b, s) in impact parameter space as 

A{s,t)=4s J d 2 b a(b,s)e iqb , (3.1) 

where q is the transverse momentum transfer vector, with, in the high energy limit, 
q 2 = —t. In this limit, a(b, s) can be assumed to be purely imaginary and therefore 
be expressed in terms of a real eikonal function x(b, s), as 

a(b, S ) = l[e-^) -1] . (3.2) 
1% 

Using ( [3.1|) and (|3.2j ) the total cross section for pp — > X can be expressed as 



a tot = -Q{A{s,t = 0)} 
s 



= 2 J d 2 b [1- ] . 

The elastic cross section is then 

a c i = 4 J d 2 b |a(b,s)| 2 

' d % h_ e -X(^s) 1 2 _ 



(3.3) 



(3.4) 



The inelastic cross section thereby reads 

Cinel = Ctot _ °"cl 



J d 2 b [1 - e~ 2 ^ ] 



(3.5) 



The elastic slope parameter at zero momentum transfer is also calculable within this 
framework and yields [20] 

B = — [ d 2 b b 2 [1 - 1 . (3.6) 

Otot J 
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3.1 Multi-parton scattering model 

The preceding expressions are completely general and model-independent, but we 
now introduce an explicit model [8-11] to predict the form of the eikonal function, 
x(h, s). First we assume that it can be decomposed into the sum of independent soft 
and hard parts, 

X tot (b, s) = XqcdO, s) + XsoftO, s) , (3.7) 

and start by considering the hard part. We consider a model in which partons are 
distributed across the face of the colliding hadrons with a spatial distribution that is 
independent of their longitudinal momentum. We assume that pairs of partons in the 
colliding hadrons scatter with independent probabilities, leading to the distribution 
of number of scatters at fixed impact parameter obeying Poisson statistics. We 
further assume that any hadron-hadron collision in which there is an elastic parton- 
parton collision above some cutoff p™ m will lead to an inelastic hadronic final state. 
Finally, we require that the inclusive cross section for hadron-hadron collisions to 
produce partons above pf in be equal to the inclusive partonic cross section folded 
with standard inclusive parton distribution functions, as given by the factorization 
theorem. The result of this model is that the inelastic cross section is given by an 
expression identical to Eq. [37|, but with x replaced by 

XQCD(b, s) = ~ A(b) 4ard(s; Pf D ) , (3-8) 

where A(h) describes the overlap distribution of the partons in impact parameter 
space and u^ rd denotes the inclusive cross section above a transverse momentum 
cutoff p t > p™ m , given by the perturbative result 

<As;pD = J2[ dx i d ^ fM fiM L ^ d&ij( ff jS \ (3.9) 

where d&ij(s)/dpt denotes the inclusive cross section for partons of types i and j and 
invariant-mass-squared s to produce jets of a given p t . 

We assume that the soft eikonal function has the same form, 

Xsoft(b, S ) = i^ ft(bK^ t , (3.10) 

where <r™£ t is the purely non-perturbative cross section below p™ m , which is a free 
parameter of the model. That is, we assume that soft scatters are the result of 
partonic interactions that are local in impact parameter. 

The elastic slope parameter discussed above relates to bulk interactions of the 
proton. Thus it can be taken as directly constraining the matter distribution "seen" 
by soft scatters. Higher p t scatters might be expected to see a different matter 
distribution, for example they might probe denser "hot spots" within the proton. 
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However, at present we take the simplest assumption for the perturbative part of the 
eikonal function, i.e. that the semi-hard scatters "see" the same matter distribution 
as the soft ones and therefore take A(b) = A so f t (h). This is clearly a simplifying 
assumption, but a reasonable one until proven otherwise. 

It is worth mentioning that according to the definition in Eq. |3.5| , the inelastic 
cross section contains all inelastic processes, including diffractive dissociation. This is 
consistent with the calculation of the inclusive hard cross section, which is calculated 
from the conventional parton distribution functions, which describe the inclusive 
distribution of partons in a hadron whatever their source, i.e. whether the proton 
remains intact or not. 



3.2 Overlap parameterization 

The only remaining freedom in the eikonal model is the functional form of the over- 
lap function A(\b\ = b). A{b) is the convolution of the individual spatial parton 
distributions of the colliding hadrons, 



A(b) 



d 2 b' G hl (\h'\)G h2 (\b-b'\). 



(3.11) 



In Refs. [9-11], G(b) is taken to be proportional to the electromagnetic form factor, 

d 2 k 



G p (b) = G p (b) 



ikb 



(2tt) 2 (l + k 2 //i 



2\2 



(3.12) 
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EM form factor, \l 2 = 2 GeV 2 
double-Gaussian, u 2 = 2 GeV 2 
EM form factor, u 2 = 0.71 GeV 2 E 
double-Gaussian, u 2 = 0.71 GeV 2 -^ 




15 



impact parameter b [\/mb] 
Figure 2: A(b) for the two parameterizations. 



\x is the only parameter and has the 
dimensions of an inverse radius. In ep 
scattering its value was measured to 
be /i 2 = 0.71 GeV 2 . This is a loose 
constraint, since the distribution of 
partons may not necessarily coincide 
with the distribution of electromag- 
netic charge. Actually, using the re- 
sults from the previous section, the 
CDF data on the total cross section 
(otot = 81.8 ± 2.3 mb [13]) and the 
elastic slope (B = 16.98 ±0.25 GeV" 2 
[19]) one can solve for the total inclu- 
sive cross section and for /z 2 , yielding 



H 2 = 0.56 ± 0.01 GeV 2 . 

In order to investigate the de- 
pendence on the assumed shape of the matter distribution, we have compared our 
default results with those obtained with a double-Gaussian distribution, as chosen 
in Refs. [21-23], 
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Gib) 



1-/3 




(3.13) 



7rr 2 



e 



7r(fc ■ r) 



Here /?, and r are all free parameters, but we choose to fix (3 and k at values that 
are reasonably generic, but also close to the tuned values used in [21-23], with the 
relative strengths given by (3 = 0.5 and the relative widths by k = 2, and view r as the 
only free parameter. The parameters /i 2 and r in the two models are arbitrary and 
should ultimately be fit to data. However, in order to have a like-for-like comparison, 
we choose to relate them in such a way that the rms value of G(b) is identical. That 
is, we describe the double- Gaussian also as being a function of /i 2 , with r set via 
brms- We illustrate the shapes of the two resulting overlap functions for two different 
values of /i 2 in Fig. 0. 

We find that for small values of fi 2 (~ 1 GeV 2 ) the results of the two models 
are extremely similar, differing at most by ±2 %. For large values (~ 3 GeV 2 ) they 
differ more, the double Gaussian distribution giving a larger de-eikonalized cross 
section (see next section) by between 30 % with the standard Donnachie-Landshoff 
total cross section prediction at the LHC and 150 % with the soft+hard Pomeron 
prediction. While these lead to somewhat different predictions, in our final results 
they effectively correspond to a distortion of the fj, 2 axis. The effect on our final plots, 
Figure |6|, is small, since our consistency requirement is mainly active at small /i 2 . 

3.3 Connection to the total cross section 

For a given point in the parameter space (p™ m , /J, 2 ) of our model, we are able to 
calculate Xqcd using Eq. [3J| The remaining uncertainty is in <j^ rd (s; p™ m ) , which 
depends on the PDF choice, the treatment of a s etc. If we now choose a value for 
the non-perturbative cross section below pf m , cr]^ t , we have the full eikonal function 
at hand and can calculate the total cross section from Eq. |3.3| . 

We will, however, turn this argument around and use the value of the total cross 
section as input to fix the additional parameter in our non-perturbative part of the 
eikonal function (o"™£ t ). For energies at which there are not yet measurements, we 
use the parameterizations of [3,4] to give an indication of the likely range of the total 
cross section. We will extract the sum o~™ rd + a^ t = Ode-eik from Eq. [O] and call 
this cross section the de-eikonalized cross section. That is, the de-eikonalized cross 
section is given by the solution to 



for a given value of /i 2 and a given value of the total cross section, a t ot- Clearly, 
Odc-cik is a function only of these two quantities. Since <J^ rd is p™ m -dependent, this 
implies that the value of a^ t we extract by this procedure is also p™ in -dependent 




(3.14) 
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v's(GeV) 

Figure 3: Cross sections in millibarn as a function of the centre-of-mass energy in GeV. 
The three different plots vary the value of fi 2 from 0.6 to 3 GeV 2 . The black curves show de- 
eikonalized total cross sections. We use the total cross section parameterization of Ref. [3] 
for the dashed curves. The solid curves use the same exponent, but the normalization 
is rescaled to fit the total cross section measurement of CDF. The dotted curve uses the 
parameterization of Ref. [4]. The coloured solid curves show cr^ ld for different values of 
pf 1111 . The coloured dash-dotted curves incorporate the simple small- x deviations discussed 
in Sect, p] 

(p™ m is a matching scale between the two sub-process cross sections and the sum of 
the two is independent of it). 

In Fig. ^| we plot the de-eikonalized cross sections for the three different total 
cross section extrapolations as a function of the centre-of-mass energy. Furthermore 
we show the value of it using different cutoffs, cr^L is now given by the difference 
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of these curves. If we interpret <7™f t as a physical cross section (the inclusive cross 
section for two partons to undergo a non-perturbative soft scattering), it cannot be 
negative. Thus the implication is that whenever the inclusive hard cross section is 
larger than the de-eikonalized one, the model is inconsistent. We will investigate this 
behaviour in more detail in Sect. PO. 



From the plots in Fig. the values for a^ t can in principle be read off. However, 
due to the logarithmic scale it is not easy to see what is implied for the energy depen- 
dence of the soft cross section. Therefore, for selected points in parameter space, a^ t 
is shown separately in Fig. |4]. Note that where the inclusive hard cross section line 
for p™ in = 3.0 GeV crosses and re-crosses the total cross-section extrapolation in the 
top left plot of Fig. |3|, the soft cross section in the top left plot of Fig. | first becomes 
negative and then positive again. The dependence of a^ t on the centre-of-mass en- 
ergy reveals two main points: First, it is noticeable that one observes a more-or-less 
constant behaviour with increasing energy only in a limited range of our parameter 
space. This behaviour is mainly present for lower values of fi 2 . Second, for the most 
extreme total cross section prediction, a^ rd is never sufficient to explain the strong 
rise with energy. There, essentially all parameter choices require a strongly rising 
soft cross section, in addition to the expected strong rise in the hard cross section. 
This is, at the very least, counter-intuitive, and one might conclude that, within our 
model, the rise of the cross section in parametrization (2.3) is too extreme. 

3.4 Parton saturation physics 

The main motivation for allowing multiparton scatterings is the rise of the inclusive 
cross section, for fixed p™ m , with increasing centre-of-mass energy. Multiparton scat- 
tering provides a mechanism to allow this quantity to exceed the total cross section. 
Eikonal models that incorporate this fact unitarize the inclusive cross section. There 
is however a second source of unitarization, the physics of parton saturation, which 
is a competing effect. To estimate the influence on our studies, we have used a simple 
modification of the PDFs recently introduced [24] into Herwig++ to mimic parton 
saturation effects for any PDF. The modification replaces xf(x) below xq by 

/ x \ Exp 

x f( x )-*\ — ] x of( x o) Vx<x , (3.15) 
where 1X0 and [Exp] are changeable parameters. HERA data indicate that saturation 



is unlikely to be a strong effect above x ~ 10~ 4 . Therefore, the strongest reasonable 
influence from this effect is obtained by setting x$ = 10 -4 , Exp = 0. The results are 
shown in Fig. |||, where the effect is visible, but small, at LHC energies. 

4. Parameter space constraints from data 

In discussing the de-eikonalized cross section, we noted that for some parameter 
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Figure 4: cr^ t for four different points in parameter space. As explained in the text, 
the extracted value of a^ t depends on the values of o" tot , fi 2 and (through the fact that 
°soft = ^de-eik — °"hard) pT 111 ■ Each panel shows a different pair of /i 2 and pf m parameters, 
while the three different curves in each use the three parameterizations for the total cross 
section as a function of energy. 

values the hard partonic cross section exceeds it. This implies in our model that the 
soft cross section should be negative. We take this as an inconsistency that would, 
for a given measured er to t at the LHC, rule out such parameter space points. In this 
section we discuss the extent to which the space of parameter values can be limited 
by this and other constraints. 

4.1 Consistency 

The parameter space in fi 2 and p™ in is shown in Fig. [5] for the Tevatron, and Fig. |G] 
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for the LHC. 

The horizontal band shows the range of /i 2 values allowed for a given value of 
the elastic slope in conjunction with the indicated cr tot . For Tevatron energies, both 
B and a to t are chosen according to the CDF measurement from Refs. [13] and [19] 
respectively. 

Our expectations on the value of the elastic slope at LHC energies simply reflect 
the range of predictions that the models of [3,4, 14, 15] give, as discussed in Sect. [2| 
For the value of ctot at the LHC we show a range of possible values motivated by 
these parameterizations. 

For a particular value of a tot (or for a given range of possible values at the LHC), 
we are able to extract constraints on the allowed parameters, by simply requiring a 
sensible performance of the eikonal model. The most basic requirement, which was 
just mentioned, is that the non-perturbative cross section that is needed to match 
the total cross-section prediction is positive. A negative value means that the model 
cannot be applied and therefore this requirement puts a stringent limit on the allowed 
values of /i 2 and pf 1111 . This limit will depend on the value of cr^d, which is not a 
stable prediction itself. We therefore calculate this limit with several variations. We 
use three different PDF sets [25-27], vary the running of a s from 1-loop, which is 
the default in Herwig+- 1- to 2-loop and finally apply the modifications to the PDF's 
described in Sect. |3.4j . The solid lines in Figs. [5] and || show these limits, where the 
entire range below the curves is excluded. The limits impose a minimal // 2 for any 
given value of p° im . Points on that line are parameter sets where <7™f t = mb. 

Another, weaker, consistency constraint we apply is related to the simulation of 
the final state of these collisions. We observe that when we embed them into the 
full simulation of Herwig++, including backward evolution of the initial state, each 
collision consumes, on average, about a tenth of the available total energy, so that 
the approximation that individual hard scatters are independent must break down, 
at least due to energy conservation, when there are more than about ten of them. 
We therefore indicate on Fig. |^ the points in parameter space where the average 
multiplicity of scatters above pf 1111 reaches 10. This is certainly not a stringent limit 
but a sensible parameter choice most likely avoids this region. 

The classic CDF analysis of the distribution in azimuth of the mean charged mul- 
tiplicity and scalar p t sum as a function of the transverse momentum of the leading 
jet [6] also provides constraints on the model as embedded in Herwig+-h Refer- 
ence [11] described the implementation of multiparton scattering into Herwig++ 
(i.e. the simulation of the final state corresponding to cr^d) and made a two- 
parameter fit (/x 2 and p™ in ) to these data. Since Herwig++ does not yet include 
a simulation of the final state corresponding to 0"™f t , we do not take the results of 
this fit as a strong constraint on the parameter space, but an indication of the effect 
such a tuning could have once a complete description is available. The result is that, 
although one obtains a best fit with the values /i 2 = 1.5 GeV 2 , p™ in = 3.4 GeV, the 
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Figure 5: The parameter space of the eikonal model at Tevatron energies. The solid 
curve imposes a minimum allowed value of /i 2 , for a given value of p™ in by requiring 
a positive value of cr"^ t . The horizontal lines correspond to the measurement of B = 
16.98 ± 0.25 GeV -2 from CDF [13]. The excluded regions are shaded. The dashed lines 
indicate the region of preferred parameter values for a fit to Tevatron final-state data from 
Ref. [11], which used the MRST2001 PDF set. The left plot shows the PDF uncertainty by 
varying the PDF set. The right plot shows the uncertainty that is implied by using 2-loop 
a s running and using the saturation modifications. 

best-fit values of the parameters are strongly correlated, with the \ 2 function having 
a long, thin, rather flat valley running from (pf 1111 = 2.5 GeV, /i 2 ~ 0.7 GeV 2 ) to 
(p™ m = 4.5 GeV, /i 2 ~ 2.5 GeV 2 ), and beyond. For any given value of pf in in this 
range one can find a /i 2 value that gives a good description of these data. 

Combining these constraints at the Tevatron, a small allowed region remains 
around pf n = 2.3 GeV and /i 2 = 0.6 GeV 2 . 

At the LHC, this region would be ruled out for all the values of er to t we have 
considered. Note that if the LHC measurement were as high a 164 mb, this would on 
its own imply an energy-dependent fj 2 , in contradiction with our initial assumptions. 

In the next section we discuss different ways in which the assumptions of the 
model might be modified to account for this potential inconsistency. 

4.2 Extensions to the model 

Some authors have suggested, within multiparton scattering models, that the pa- 
rameters of the model, analogous to our /x 2 and pf" n , should be energy dependent. 
In this section we briefly discuss the arguments for these models. 

In [28] a simple model of the spatial/momentum structure of a hadron was con- 
structed from which it was argued that the colour screening length decreases slowly 
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with increasing energy. This translates into a p™ m that increases slowly with energy, 
for which they estimated p™ m ~ s e with e in the range 0.025 to 0.08. The actual value 
used in Refs. [21-23] is 0.08, leading to a 35% increase in p™ m from the Tevatron to 
the LHC. 

In [29,30] a multiparton model was constructed that is very similar to ours at 
low energy, with an impact parameter distribution of partons given by the electro- 
magnetic form factor. However, soft gluon effects were estimated and summed to all 
orders, to give a mean parton-parton separation, b rms , that falls with energy, quickly 
at first, but then saturating: the value at 1 TeV is about a factor of two smaller than 
at low energy, while the value at 14 TeV is only about 10% smaller still. In terms of 
our simple model in which the matter distribution always has the form factor form 
and is parameterized by fx 2 , (b 2 ) oc l//i 2 and this corresponds to jj 2 ~ 2.8 GeV 2 at 
the Tevatron and ~ 3.4 GeV 2 at the LHC. Not only would this introduce an energy 
dependence in /i 2 , but the values imply a different /i 2 for hard partonic interactions 
than that derived from the measured elastic slope parameter, a point that we will 
return to below. 

Note that both these sources of energy dependence would act in the right di- 
rection to evade the potential consistency constraints at the LHC. Allowing p™ m 
and/or /i 2 to increase with energy would move the model towards the allowed region 
in Fig. ||. 

The CDF collaboration have published measurements of the double-parton scat- 
tering cross section [31,32]. As pointed out in Ref. [33] the quantity called cx eff there 
is not the effective cross section as it is usually defined, 

0"eff = o, (4.1) 

fd 2 h(A(b)) 2 ' 

but is related to the latter by a small correction. Using the value of this correction 
estimated in Ref. [34], we obtain <j c g ~ 11.5 ± 2 mb. In our form factor model, this 
corresponds to /i 2 ~ 3.0 ± 0.5 GeV 2 . It is interesting to note that this value is close 
to the one predicted by the analysis of Refs. [29,30] mentioned earlier. Again, this 
value is inconsistent with the assumption that the hard scatters "see" a form factor 
matter distribution derived from the elastic slope parameter. Recall from our earlier 
discussion that we do not expect significant qualitative differences for other models of 
the matter distribution in the proton, merely some distortions of the parameter-space 
plane. 

Improved analyses of these and other observables are under way and, once com- 
pleted, in particular with a simulation of the final state of cr™f t , will provide strong 
constraints on the values of the parameters fx 2 and p™ in in our model. 
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5. Conclusions 

The connections between our underlying event model and the total proton-proton 
cross-section have been discussed. Requiring consistency of the model up to LHC 
energies imposes constraints on the allowed parameter values, for a given range of 
possible measurements of otot at the LHC. Our main result is summarized in Fig. || 
which shows these constraints for various values of the total cross section at the LHC 
and various inputs to the perturbative cross section calculation. Taking the Tevatron 
data together with the wide range of possible values of cr tot considered at LHC, no 
allowed set of parameters (/i 2 , pf m ) remains for our simple model. 

This would imply that soft and hard scatters see different matter distributions as 
a function of impact parameter and/or that the parameters of our model are energy 
dependent. The measurement of double-parton scattering at the Tevatron supports 
the idea that hard scatters see a more dense matter distribution than is implied by 
the t-slope of the elastic cross section. Various phenomenological models also predict 
such effects. 

With steadily improving data from the Tevatron, more detailed phenomenologi- 
cal analyses being performed and the prospect of data from the LHC, there is a real 
hope that the various simplifying assumptions that go into our model will be tested to 
the limit and we will discover where, if anywhere, more detailed understanding of the 
dynamics of underlying event physics is needed. The biggest unknown in our analysis 
is the total cross section at the LHC. With even a first imprecise measurement of 
this cross section, we could strengthen our parameter space analysis considerably. 
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Figure 6: The parameter space of the eikonal model and three constraints. The first one 
drawn as solid curve imposes a minimum allowed value of /i 2 , for a given value of pf 1111 
by requiring a positive value of c^oft- The second one, in dashed lines is deduced from an 
upper limit of the average number of additional semi-hard scatters in a typical minimum 
bias event. The excluded regions are shaded. The third constraint comes from the expected 
range of values for the elastic slope B. From top to bottom we calculate these constraints 
for the three different total cross sections at LHC, discussed in Sect. [2|, always with the 
same range of B = 19 — 22 GeV -2 . Finally the left column shows the PDF uncertainty 
by varying the PDF set. The right column shows the uncertainty that is implied by using 
2-loop a s running and using the saturation modifications. 



